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ABSTRACT 
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, Aims. This paper reports the photospheric, magnetic and circumstellar gas characteristics of the magnetic B8p star 36 Lyncis 
23 ■ (HD 79158). 

Methods. Using archival data and new polarised and unpolarised high-resolution spectra, we redetermine the basic physical 
<*■ . properties, the rotational period and the geometry of the magnetic field, and the photospheric abundances of various elements. 
J s J~j ' Results. Based on magnetic and spectroscopic measurements, we infer an improved rotational period of 3.83475 ± 0.00002 d. We 
determine a current epoch of the longitudinal magnetic field positive cxtrcmum (HJD 2452246.033), and provide constraints 
■ on the geometry of the dipole magnetic field (i > 56°, 3210 G < Bd < 3930 G, /3 unconstrained). We redetermine the effective 
temperature and surface gravity using the optical and UV energy distributions, optical photometry and Balmer line profiles 
(T e ff = 13300 ± 300 K, \ogg = 3.7 — 4.2), and based on the Hipparcos parallax we redetermine the luminosity, mass, radius 
and true rotational speed (L = 2.54 ± 0.16 L Q ,M = 4.0 ± 0.2 M Q ,R = 3.4 ± 0.7 R Q ,v cq = 45 - 61.5 kms^ 1 ). We measure 
photospheric abundances for 21 elements using optical and UV spectra, and constrain the presence of vertical stratification 
of these elements. We perform preliminary preliminary Doppler Imaging of the surface distribution of Fe, finding that Fe is 
distributed in a patchy belt near the rotational equator. Most remarkably, we confirm strong variations of the Ha line core 
which we interpret as due to occultations of the star by magnetically-confined circumstellar gas. 



Key words, stars: individual: 36 Lyn - stars: magnetic fields - polarisation 



1. Introduction 
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* Based on observations obtained using the MuSiCoS spec- 
tropolarimeter at Pic du Midi Observatory, the Canada- 
France-Hawaii Telescope (operated by the National Research 
Council of Canada, the Centre National de la Recherche 
Scientifique of France, and the University of Hawaii), the 
Dominion Astrophysical Observatory, Herzberg Institute of 
Astrophysics, National Research Council of Canada, the David 
Dunlap Observatory, University of Toronto, and the Elginfield 
Observatory, University of Western Ontario. 



The magnetic Bp stars exhibit a rich variety of magnetic 
and wind phenomena, many of which are unique to this 
class of star. Landstreet & Borra (1978), detecting for the 
first time a magnetic field in a helium strong star, pro- 
posed that the photometric and spectroscopic variations 
of a Ori E resulted from hot gas trapped above the mag- 
netic equator. Brown et al. (1985) discovered unexpect- 
edly strong C IV, Si iv and N v ultraviolet resonance 
line absorption and emission that varied with the rota- 
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tional period of the magnetic helium weak star HD 21699. 
Subsequent work revealed similar phenomena in a num- 
ber of other Bp stars, which was attributed to rotational 
modulation of magnetically-structured stellar winds (e.g. 
Shore 1987, Shore et al. 1990, Shore & Brown 1990, Bolton 
1994, Donati et al. 2001, Smith & Groote 2001, Neiner et 
al. 2003). The additional discovery of non-thermal radio 
emission (Drake et al. 1987) from a subset of these objects 
resulted in the proposal that these stars are surrounded 
by an extended wind-generated magnetosphere (Linsky et 
al. 1992), in which the radio luminosity results from opti- 
cally thick gyrosynchrotron emission and should therefore 
correlate with both magnetic field strength and effective 
temperature. According to Drake et al. (2002), by 2002 
about 120 magnetic CP stars had been observed in the 
radio, and of these about 25% had been confidently de- 
tected. The general existence of this phenomenon, and the 
establishment of correlations between effective tempera- 
ture, magnetic field strength, rotational period and radio 
emission suggest that this phenomenon is intrinsic to Ap 
and Bp stars, and is related to their magnetic properties. 
Subsequently, X-ray emission was discovered (Drake et al. 
1994) from a small number of these objects. As pointed 
out by Drake et al. (2002), no clear correlation of X-ray 
brightness with other physical or observational character- 
istics (including radio luminosity) has been established, 
suggesting (at least in the A-type magnetic stars) that 
X-ray emission may be due to the presence of late-type 
companions, and unrelated to their magnetic properties 1 . 
More recently, sophisticated models describing the hydro- 
dynamic interaction between the magnetic field and wind 
flow have been developed (Babel & Montmerle 1997, ud- 
Doula & Owocki 2002, Preuss et al. 2004, Townscnd & 
Owocki 2005), leading to the possibility of combining the 
various individual kinds of observations (magnetic, radio, 
optical, UV, X-ray) into a single coherent picture of the 
wind-magnetic field interaction and subsequent circum- 
stellar structure of magnetic Bp stars. 

36 Lyncis (HD 79158) is a bright (V = 5.3), well- 
studied magnetic (Borra, Landstreet & Thompson 1983) 
helium weak B star. It is classified as B8IIIpMn in the 
Bright Star Catalogue (HofHeit et al. 1995), as a Si star 
by Molnar (1972), as a Mn star by Cowley (1972), and as 
a Sr Ti helium weak star by Borra & Landstreet (1983). 
It was first identified as peculiar by Edwards (1932). The 
spectrum and abundances have been studied by Searle & 
Sargent (1964), Mihalas & Henshaw (1966) and Sargent, 
Greenstein & Sargent (1969). According to these studies, 
lines of neutral He and O are weak, whereas those of C, 
Ne, Si, Mg and Cr are only marginally anomalous. Lines 



1 There exists compelling evidence, in the form of variability 
of the X-ray flux according to the stellar rotational period, that 
some magnetic B and O stars (the Bp star AB Aur: Babel & 
Montmerle (1997); the magnetic 07V star 9 1 Ori C: Gagne et 
al. (1997), Donati et al. 2002) are intrinsic X-ray emitters. For 
the majority of magnetic upper-main sequence stars, however, 
this is not the case. 



of Ti and Fe are anomalously strong. HD 79158 is a mem- 
ber of the sn class (Abt & Levato 1978), which means 
its Si II lines are sharp, whereas He I lines are nebulous. 
Ryabchikova & Stateva (1995) concluded that 36 Lyn is 
a classical helium weak star. 

Using IUE spectra, Sadakane (1984) reported the pres- 
ence of strong lines coincident with the position of the 
C iv AA1548, 1550 and Si iv AA1394, 1403 resonance dou- 
blets. Shore, Brown & Sonneborn (1987) obtained addi- 
tional magnetic field measurements and IUE spectra of 
36 Lyn and two additional sn stars, confirming the pres- 
ence of the "superionised" C iv and Si iv resonance lines 
and detecting variability of these features. Shore et al. 
(1990) extended these observations and analysis, obtain- 
ing 22 new measurements of the longitudinal magnetic 
field and determining a rotational period of 3.8345 ±0.001 
days. 

The goal of the present paper is to provide an accurate 
reassessment of the magnetic and photospheric properties 
of 36 Lyn as the basis for a detailed study of properties of 
the circumstellar matter (Smith et al. 2006, "Paper II"). 

2. New and archival observations 

2.1. Un polarised optical spectra 

The journal of the unpolarised optical spectroscopic obser- 
vations is presented in Table 1 , available only in electronic 
form from the CDS. 

2.1.1. DAO and DDO Spectra 

Ha observations of 36 Lyn were obtained at the Dominion 
Astrophysical Observatory (DAO) with both the 1.85 m 
Plaskett (1991, 1992, and 1997-2000) and 1.22 m McKellar 
(2001, 2003-2005) telescopes, and at the David Dunlap 
Observatory 1.88 m telescope (1992, 1997, and 2004). 

At DAO, on the 1.85m we used the Cassegrain spec- 
trograph with the 21 inch //5 camera, 1800 line mm -1 
grating in first order, VSIS21B image sheer, and various 
CCDs. The resulting spectra covered a minimum of ap- 
proximately 100 A and, when possible, included the Hei 
A6678 feature. The dispersion of 10 A mm -1 and effec- 
tive slit width of 40 /im gave a resolution of approximately 
R = 16400. 

The 2001 McKellar observations were obtained 
on the coude spectrograph with the short camera, 
1200 line mm -1 holographic grating in first order, IS32R 
image sheer, and the UBC-1 CCD. These short camera 
observations have a dispersion of 10.1 Amm -1 and the 
resolution is approximately R = 20000. 

For the 2003, 2004 and 2005 McKellar data we again 
used the coude spectrograph, but with the long camera, 
830 line mm -1 grating in first order, IS32R image sheer, 
and the SITe4 CCD with 15 (im pixels. The dispersion of 
4.8 Amm -1 and projected slit width of 36 /xm produced 
a resolution of R = 38000. 
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Exposure times on both telescopes ranged from 10 to 
30 minutes. 

DDO spectra obtained in 1992 and 1997 were recorded 
with a Thomson 1024x1024 CCD on the Cassegrain spec- 
trograph of the 1.88 m telescope. The effective resolution 
was R = 11600. Most of the spectra were centered near 
6600 A so that the 200 A range covered by the detector 
included both Ha and the He I A6678 A line. 

The DDO spectra obtained in 2004 were recorded with 
a Jobin-Yvon 800x2000 thinned, back illuminated CCD 
with 15 fim pixels. The effective resolution of these spectra 
is 14700. These observations were taken in the same way 
as the earlier observations. However, there were difficulties 
because the detector system was not operating properly. 
Due to the large RFI introduced by the heater system 
and the lack of a cold finger in the dewar, we had to op- 
erate the detector system with the temperature floating 
near the liquid nitrogen vaporization temperature. As a 
result, there were large changes in the CCD temperature 
as a function of hour angle of the telescope. To minimize 
this problem, we took frequent biases, and linearly inter- 
polated in time between them to obtain the bias for each 
stellar spectrum. The low operating temperature also re- 
duced the dynamic range of the CCD. Flat field exposures 
showed that the response started to become nonlinear at 
24k e~. As a result, these spectra have somewhat lower 
S/N than our earlier spectra. 

All data were reduced with IRAF 2 in a conventional 
manner, including bias subtraction, flat-fielding, wave- 
length calibration, continuum normalisation and removal 
of telluric absorption lines. 
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Fig. 1. LSD mean Stokes IV QU profiles. Rotational 
phases (according to the ephemeris given in Eq. (2) and 
multiplied by 1000) are given in the left panel. Apart from 
a marginal detection at phase 0.897, Stokes Q and U pro- 
files are consistent with zero linear polarisation across the 
mean line. Variability of the Stokes I profile is observed, 
and modulation of the Stokes V profile shape by an ap- 
parently non-uniform Fe surface abundance distribution 
(analogous to that encountered for Aur by Wade et al. 
2000a, for example) supports the view that variability of 
the mean Fe line results primarily from non-uniform sur- 
face abundances (see Sect. 6). 



2.1.2. CFHT Spectra 

Canada-France-Hawaii Telescope (CFHT) data were ob- 
tained on 1991 November 19-22 and 1995 October 3 with 
the coude spectrograph and red mirror train, //8.2 spec- 
trograph, 600 line mm -1 grating in first order, and red 
Richardson image sheer. Detectors for these runs were the 
1872 diode Reticon array (1991) and Loral3 (1995) CCD. 
The effective resolution was approximately R = 22000. 
Exposure times were approximately 15 minutes with the 
Reticon and 10 minutes with the CCD. A thorium-neon 
hollow cathode lamp was used as a comparison line source 
for the Reticon data while a thorium-argon lamp was used 
for the CCD observation. In the case of the Reticon data, 
each spectrum (including stellar, comparison lamp, and 
flat-field exposures) consisted of a raw exposure and eight 
1 second baseline exposures and the first step of the data 
reduction was the subtraction of the average of these base- 
lines from each raw spectrum. Flat-fielding, wavelength 
calibration, continuum rectification, and telluric line re- 
moval was then carried out as described for the DAO ob- 

2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., under 
cooperative agreement with the National Science Foundation. 



servations above. The CFHT CCD observations were also 
reduced in an identical fashion to the DAO data. 

Heliocentric velocity corrections have been applied to 
all of the unpolarised spectra. 

2.2. Circular and linear polarisation spectra 

Circular polarisation (Stokes V) and linear polarisation 
(Stokes Q and U) spectra of 36 Lyn were obtained dur- 
ing Feb/Mar 2000 and Dec 2001 using the MuSiCoS spec- 
tropolarimeter mounted on the 2 metre Bernard Lyot tele- 
scope at Pic du Midi observatory. The spectrograph and 
polarimeter module are described in detail by Baudrand 
& Bohm (1992) and by Donati et al. (1999) respectively. 
The standard instrumental configuration allows for the 
acquisition of circular or linear polarisation spectra with 
a resolving power of about 35,000 throughout the range 
4500-6600 A. 

A complete Stokes V circular polarisation observation 
consists of a series of 4 subexposures between which the 
polarimeter quarter-wave plate is rotated back and forth 
between position angles of —45° and +45°. This procedure 
results in exchanging the orthogonally polarised beams 
throughout the entire instrument, with the aim of re- 
ducing systematic errors in polarisation measurements of 
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Table 2. Journal of MuSiCoS spectropolarimetric observations. Phases correspond to the midpoint of the Stokes V 
exposure, are calculated according to the new ephemeris described in the text, and exposure times are in seconds. 
Reported S/N are peak S/N per 4.5 kms -1 pixel. TVlsd is the noise level in the associated LSD profile, per 4 kms -1 
LSD velocity bin. 



sharp-lined stars to below a level of about 0.01% (Wade 
et al. 2000a). The procedure for recording Stokes Q and U 
linear polarisation observations is similar, except that the 
entire instrument is rotated (rather than the quarter-wave 
plate, which is removed from the beam for linear polari- 
sation observations). Further details about the observing 
procedure are provided by Wade et al. (2000a). 

In total 11 Stokes V, 8 Stokes Q and 8 Stokes U spec- 
tra of 36 Lyn were obtained, with peak signal-to-noise 
ratios (S/N) of typically 300:1 per pixel in Stokes V and 
325:1 in Stokes Q/U. In addition, during each run obser- 
vations of various magnetic and non-magnetic standard 
stars were obtained within the context of other observing 
programmes (e.g. Shorlin et al. 2002, Kochukhov et al. 
2004) which serve to verify the nominal operation of the 
instrument. 

The log of spectropolarimetric observations of 36 Lyn 
is shown in Table 2. 

2.3. IUE spectra 

The IUE spectra employed in this study were the 
NEWSIPS-extracted, large aperture high-dispersion files 
from the Short Wavelength Prime (SWP) camera. These 
were downloaded from the MAST data archive. The data 
used were all 24 IUE SWP/large- aperture spectra for 36 
Lyn, obtained during 3 epochs in 1985, 1987 and 1988. 
The data are described in more detail by Smith & Groote 
(2001). 

2.4. Balmer-line Zeeman analyser polarimetric 
observations 

Three measurements of the longitudinal magnetic field 
of 36 Lyn were obtained in 1992 using the University 
of Western Ontario (UWO) photoelectric polarimeter 
mounted on the 1.2 m telescope at the UWO Elginfield ob- 
servatory. The polarimeter was employed as a Balmer-line 



Zeeman analyser to measure the fractional circular polar- 
isation in the wings of H/3 at ±5.0 A from line centre. 
The longitudinal magnetic field was inferred from mea- 
sured fractional circularly polarised flux using a conver- 
sion factor of 17500 G/%, as determined from scans of the 
H/3 profile of 36 Lyn. The instrument and observing tech- 
nique are described in detail by Landstrcet (1982). The 
measurements are reported in Table 3. 

3. Least-Squares Deconvolution and longitudinal 
magnetic field 

We begin by applying the Least-Squares Deconvolution 
(LSD) multiline analysis procedure (Donati et al. 1997) to 
each of the MuSiCoS polarised spectra in order to obtain 
mean Stokes IQUV profile sets (e.g. Wade et al. 2000a). 
LSD mean profiles were extracted using line masks derived 
from VALD line lists, corresponding to T e s = 13500 K, 
\ogg — 4.0 (see Sect. 5), a minimum line depth of 10% and 
an abundance table specific to 36 Lyn (using abundances 
derived in Sect. 6). More information about the extraction 
of LSD profiles and the construction of LSD line masks is 
provided by Wade et al. (2000b) and Shorlin et al. (2002). 

The LSD profiles (extracted for lines of Fe) are shown 
in Fig. 1. The Stokes V profiles show complex structure, 
are clearly detected at all phases, and are strongly modu- 
lated as the star rotates. Remarkably, Stokes Q and U are 
never confidently detected, notwithstanding a typical LSD 
noise level below 0.015%. The Stokes I profiles also show 
variability, which we interpret as rotational modulation of 
a non-uniform surface abundance distribution of Fe. This 
will be discussed in further detail in Sect. 7. 

The longitudinal magnetic field (B z ) and its formal 
uncertainty gb were inferred from each of the extracted 
LSD Stokes I/V profile sets by numerical integration over 
velocity, in the manner described by Wade et al. (2000b). 

The 11 new LSD longitudinal field measurements have 
typical uncertainties <jb ~ 60 G, and range from -790 G 
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2448739.621 


0.622 


-90 ± 200 


2448744.609 


0.922 


-790 ± 320 


2448754.676 


0.548 


660 ± 260 


2451578.555 


0.940 


-612 ± 48 


2451584.425 


0.470 


595 ± 66 


9451 586 556 


0.026 


—91 ± 109 


2451596.520 


0.624 


24 ± 56 


2451597.766 


0.949 


-762 ± 52 


2451601.519 


0.928 


-694 ± 68 


2451603.512 


0.448 


662 ± 82 


2451606.479 


0.222 


876 ± 68 


2451608.496 


0.747 


-673 ± 50 


2452245.637 


0.897 


-790 ± 90 


2452246.611 


0.151 


609 ± 87 



Table 3. New longitudinal magnetic field measurements 
of 36 Lyn obtained from Balmer-line photopolarimetry 
(HJD 2448739-2448754) and from MuSiCoS LSD mean 
Stokes I and V profiles (HJD 2451578-). Phases are cal- 
culated according to the rotational ephemeris, Eq. (1), 
determined in Sect. 4. All uncertainties are la. 



to +876 G. The inferred values of the longitudinal field of 
36 Lyn are reported in Table 3. 

We have combined the 14 new measurements of 
the longitudinal magnetic field with 24 previously pub- 
lished measurements (2 reported by Borra, Landstreet 
& Thompson (1983) and 22 reported by Shore et al. 
(1990), all obtained using a Balmer line Zeeman anal- 
yser at H/3) in order to determine the magnetic period. 
The periodogram of the longitudinal field data, obtained 
using a modifed Lomb-Scargle technique, is characterised 
by a strong, unique peak at 3.83495 ± 0.00003 days (re- 
sulting in a reduced \ 2 of 2.07 for a first-order harmonic 
fit, and where the formal la error bar is derived using 
the reduced x 2 statistical tables of Bevington 1969). This 
period is consistent with, although more precise than, 
published magnetic and photometric values of the period 
(e.g. Shore et al. 1990, confirmed by Adelman (2000), 
P rot = 3.8345 ± 0.001 days). This period has been kindly 
confirmed using minimum-false alarm probability meth- 
ods by P. Reegen, who finds a best-fit period of 3.83495 
days with a significance \ogS = 3.34. 

It should be noted that a second-order fit to the 
magnetic data provides a significantly lower reduced \ 2 
(1.31). Independent fitting of the individual datasets (the 
Balmer-line measurements and the LSD metallic-line mea- 
surements) demonstrates that the Balmer-line data are 
well-described by a first-order fit, and that the second- 
order contribution comes entirely from the LSD data. The 
shapes of the longitudinal field variation as defined by the 
two datasets are therefore quite clearly different. This is 
not too surprising - it is well documented that longitudinal 
field variations derived from metallic-line measurements 
often contain non-sinusoidal contributions (e.g. Wade et 
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Fig. 2. Variation of the Ha profile core of 36 Lyn due to 
the magnetically-confined circumstellar disk (see Sect. 4 
and 7). 

al. 2000b), likely due to the presence of non-uniform dis- 
tributions of these elements across the stellar surface. 

4. Variations of the Ha profile 

Upon examination of the optical spectra, it was immedi- 
ately noted that the core depth of the Ha line of 36 Lyn is 
strongly variable. Variability was first reported by Takada- 
Hidai & Aikman (1989), who attributed it to occultation 
of the stellar disc by the magnetically-confined circum- 
stellar plasma. These authors obtained several observa- 
tions of Ha, and noted that the core occasionally weakens, 
possibly accompanied by a depression of the inner wings. 
We confirm the variability, which is clearly intrinsic to 
Ha, and not (for example) the result of blending with 
strongly-variable metallic lines (the metallic-line variabil- 
ity of 36 Lyn is too weak, and no suitably-located lines 
exist in any case) or due to continuum registration errors 
(the continuum near Ha in our spectra is uniformly excel- 
lent, and shows internal and external agreement amongst 
the datasets to within about 1%). However, in contrast to 
the observations of Takada-Hidai & Aikman, we observe 
no variability of the inner wings, and we observe a core 
which occasionally strengthens. 

We have measured equivalent widths of the core of Ha 
in all of our spectra, in the range 6561.0 to 6565.5 A, 
after local reregistering of the spectra to further re- 
duce systematic differences in continuum normalisation 3 . 
Uncertainties associated with the equivalent widths were 
also calculated by adopting the inverse continuum signal- 
to-noise ratio as an approximate standard error bar asso- 
ciated with each spectrum pixel, and propagating these 
errors through the equivalent width calculation. We find 
that the measured Ha equivalent width varies rather 
rapidly as a function of phase, with significant variations 
of about 15% (see Fig. 2). 

When the Ha equivalent width measurements are 
folded according to the magnetic period determined in 

3 The Ha equivalent width measurements are reported in 
Table 1, available only in electronic form from the CDS. 
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Phase (P=3.83475 days) 

Fig. 3. Upper frame - Variation of the Ha profile core 
equivalent width and longitudinal field of 36 Lyn. Detail 
of primary occultation, measurements from all spectra, 
phased according to the adopted magnetic rotational pe- 
riod of 3.83495 days. Note the clear and systematic offset 
of measurements obtained at different epochs (represented 
by different symbols). Lower frame - Longitudinal mag- 
netic field measurements phased according to the equiva- 
lent width period 3.83475 days. Note the apparent offset 
between the older (H/3, open squares) and more recent 
(MuSiCoS, filled circles) longitudinal field measurements. 

Sect. 3, the various epochs are found to be offset in phase 
relative to one another by up to 0.06 cycles, as shown in 
Fig. 3 (upper frame). If we interpret these offsets to be the 
result of an error in the period determined from the mag- 
netic measurements, they imply a best-fit equivalent width 
period (determined using Phase Dispersion Minimisation; 
Stellingwerf 1978) of 3.83475±0.00002 days. When phased 
according to this period, the equivalent width data show 
a clear, coherent variation, with a sharp maximum near 
phase 0.3, as well as a much weaker increase about one-half 
cycle later (Fig. 4). We interpret these variations in the 
context described by Shore, Brown & Sonneborn (1987). 
In this context, the increase in Ha absorption results from 
occultation of the stellar disc by plasma which is con- 
fined magnetically near the magnetic equator, and forced 
to co-rotate with the star. The (neutral) H gas is coupled 



to the ionised g clS elS cl result of frequent reionisations, 
as well as by collisions. This interpretation is supported 
by the presence of two peaks, separated in phase by ap- 
proximately 0.5 cycles, and approximately coincident with 
magnetic crossover phases (when (B z ) ~ 0). It is further 
supported by observed systematic changes in the radial 
velocity (RV) of the Ha core coincident with the occulta- 
tion phases. As would be expected, the core RV is initially 
blueshifted (with respect to the systemic RV) at the be- 
ginning of the occultation, and becomes redshifted as the 
occultation ends (this is illustrated in Fig. 5). 

The intensities of the two occultation peaks are clearly 
unequal. This implies that the column depth of neutral H 
occulting the disc differs substantially at phases 0.0/0.5, 
i.e. that the occulting material is not distributed uniformly 
within or near the magnetic equatorial plane. The accumu- 
lation of gas overdensity at the intersections of the mag- 
netic and rotational equators (i.e. at the nodal points) of 
magnetic Bp stars, and differences in the column density 
of these "clouds" , has been documented for the He-strong 
star a Ori E (e.g. Landstreet & Borra 1979), for example. 

The physical conditions, structure and dynamics of the 
magnetically-confined wind will be discussed in more de- 
tail in Paper II. 

The period deduced from the Ha equivalent width 
data is not formally consistent with that determined using 
the magnetic data. However, when the magnetic data are 
phased according to the Ha period, an acceptable phase 
variation is achieved. Although the new magnetic data 
show a small apparent offset in phase with respect to 
the archival data (Fig. 3, lower frame), this offset could 
easily result from well-known systematic differences in 
the shapes of field variations diagnosed from metallic vs. 
Balmer lines. 

We therefore adopt the Ha period as the rotational 
period of 36 Lyn. Hereinafter, all data have been phased 
according to the ephemeris: 

JD = (2443000.451 ± 0.03) + (3.83475 ± 0.00002) • E, (1) 

where the zero-point corresponds to the epoch of maxi- 
mum Ha core absorption. 

5. Fundamental characteristics of 36 Lyn 

5.1. Fundamental parameters and evolutionary status 

We begin by redetermining the effective temperature and 
surface gravity of 36 Lyn using published Geneva and 
uvby/3 photometry, the observed spectral energy distribu- 
tions at optical and UV wavelengths, and the Ha and H/3 
profiles in our spectroscopic data. 

Geneva photometry (from the Geneva photometric 
database, http : / / obswww . unige . ch/gcpd/ phl3 . html 
Burki 2002) coupled with the calibration of North & 
Nicolet (1990) give T cff = 13650 K and logg = 3.75. 
Stromgren photometry (Hauck & Mermilliod 1998) cou- 
pled with the calibration of Moon & Dworetsky (1985) 
give T off = 13740 and log 5 = 3.65. 
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Rotational phase 

Fig. 4. Variation of the Ha profile core equivalent width of 36 Lyn, with measurements phased according to the 
equivalent width period 3.83475 d. Note the significant difference in amplitude of the </> ~ 0.0 and <f> ~ 0.5 peaks. Error- 
bars shown correspond to 2a confidence. 




Rotational phase 

Fig. 5. Variation of the equivalent width and radial veloc- 
ity of the Ha profile core of 36 Lyn during the primary 
occultation at phase 0.0. Measurements have been phased 
according to the equivalent width period 3.83475 d. Error 
bars shown correspond to 2a confidence. 



Using theoretical flux spectra, we have fit the observed 
optical (Adelman & Pyper 1983) and UV (Jamar et al. 
1976) spectrophotometry of 36 Lyn. Using a modified 
version of the ATLAS9 model atmosphere code (Kurucz 
1993) and our photosphcric abundance determinations 
(Sect. 6), we have performed this calculation for a grid 
of temperatures and gravities bracketing the photomet- 
ric values, and various metallicities ranging from 1 — 10 x 
solar. We have offset magnetic broadening in the atmo- 
sphere model by assuming a 2 kms -1 microturbulence. 



In Fig. 6 we show the observed energy distribution along 
with selected models. The best model fit to the energy 
distribution of 36 Lyn (T cff = 13300 K, logg = 4.0, 
[M]=0.5) is shown in Fig. 6 (upper panel), as well as 
fits for (T cff = 13600 K, logg = 3.7, [M]=0.0) and 
(T cff = 13000 K, logg = 4.0, [M]=1.0). For the UV fluxes, 
we used TD1 data which was converted from flux into 
monochromatic magnitudes. We then reduced to 5000 A 
using m(5000 A)=5.37 (a value very close to the known 
my=5.32). 

Finally, using ATLAS9 model Balmer line profiles for 
solar and enhanced metallicities, we have fit the observed 
wings of the Ha and H/3 profiles and obtain a best-fit for 
T e g = 13600 K and logg = 3.7 (Fig. 6, lower frames). 

All three models provide a satisfactory fit to the ob- 
served flux distribution in the optical region, while they 
differ in the UV region. As E(B - V) = 0.003 for 36 Lyn 
(see, for example, Lucke IT9"TBf ■ this cannot be explained 
as due to unaccounted-for reddening. The model with the 
highest metallicity, 13000g40pl0k2, provides the better fit 
belowA=1800A (1/A( microns)=5.5), while the moderate 
metallicity model, 13300g40p05k2, fits better in all spec- 
tral regions longer than 1800 A. The higher temperature 
model does not fit the UV flux distribution, but it pro- 
vides the best fit to the observed hydrogen line profiles. 
These systematic and significant discrepancies probably 
show true limitations of ATLAS9 model atmospheres in 
describing the photospheric structure of 36 Lyn. 

As a compromise, we have adopted T e g = 13300 ± 
300 K for 36 Lyn, and have chosen a model with 
T e ff =13300 K, log <7=4.0, and a logarithmic metallicity 
with respect to the sun of 0.5 for our abundance analy- 
sis. These values have been selected giving approximately 
equal weight to the energy distribution and Balmer lines, 
although we point out that they are rather arbitrary - it 
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Mass (Mo) 
Radius (72.0 ) 
Age (y) 
T eS (K) 
logs' 
v sini 
Period 
sin i 



4.0 ±0.2 
3.4 ±0.7 
0.79 - 1.1 x 10 s 
13300 ± 300 
3.7 - 4.2 
48 ±3 

3.83475 ± 0.00002 d 
0.83-1.00 



Table 4. Fundamental characteristics derived for 36 Lyn 



is really impossible to confidently select any model within 
the range 13000 < T cff < 13600 K, 3.7 < log 5 < 4.0 or 
0.0 <[M]< 1.0. 

Now, using the derived temperature, along with the 
apparent visual magnitude (my = 5.3), absolute magni- 
tude (M v = —0.8) and Hipparcos distance (d = 162 ± 
15 pc) reported by Gomez et al. (1998), we obtain the 
bolometric magnitude of 36 Lyn Mboi = — 1.65 ± 0.4 (us- 
ing a bolometric correction BC = —0.85 determined using 
the interpolation formula of Balona 1994). The luminosity 
is therefore logL/L© = 2.54±0.16, and the stellar radius 
R = 3.4 ± 0.7 Rq. From the logT - \ogL/L Q HR dia- 
gram position of 36 Lyn (illustrated in Fig. 7), we obtain 
the stellar mass M = 4.0 ± 0.2Af© and age logi = 8.0 
(t = 79 — 110 Myr) based on the model evolutionary cal- 
culations for solar metallicity of Schaller et al. (1992). The 
combination of R and M give a logarithmic surface grav- 
ity of log g = 4.0 ± 0.2. We point out, as elaborated by 
Bagnulo et al. (2005), that the use of a bolometric cor- 
rections derived for normal stars, coupled with compar- 
isons with solar metallicity evolutionary tracks, and un- 
certainties regarding the application of the Lutz-Kelkcr 
correction, likely underestimates the luminosity and age 
uncertainties. 

The measured projected rotational velocity vsini = 
48 ± 3 kms -1 , obtained using individual metallic line 
profiles (see Sect. 6), and rotational period of 3.83475 ± 
0.00002 days imply a projected radius (assuming rigid ro- 
tation) of R = (3.64 ± 0.23) sin i R Q . The radius deter- 
mined from the HR diagram is numerically smaller than 
the projected radius determined from rigid rotation (al- 
though they are consistent within the uncertainties). The 
combination of the two radius estimates implies i > 56°, 
with a best-fit value of i — 90° , and a true rotational speed 
in the range 45-61.5 kms -1 . 

The derived fundamental and rotational characteris- 
tics of 36 Lyn are summarised in Table 4. 

5.2. Binarity 

There is no indication of binarity in the literature, and the 
Hipparcos mission classifies 36 Lyn as single. The Stokes / 
LSD profiles are constant in mean radial velocity at ±22 ± 
2 kms -1 . The Ha profiles, outside of occultation phases, 
give ±23 ± 1 kms" 1 . Stickland & Weatherby (1984) find 
v la d — ±26 ± 1 kms -1 , whereas the Bright Star Catalogue 
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Fig. 7. HR diagram position of 36 Lyn. The filled circle 
represents the derived position of 36 Lyn, solid curves are 
evolutionary tracks for 3.5, 4.0 and 4.5 M Q for solar metal- 
licity (Schaller et al. 1992), and the dashed curve is the 
associated isochrone for 10 8 y. 

(HofHeit et al. 1995) gives ±21 kms" 1 , Aikman (1976) 
gives ±21.9 kms" 1 , Abt (1970) gives ±26.4, and Takada- 
Hidai & Aikman (1989; measured from Ha) give ±29.7 ± 
1.6 kms -1 . Although there exists some diversity among 
these measurements, this is expected given the metallic- 
line RV variations (due to surface features - see Sect. 6.3) 
and the Ha RV variations (due to the circumstcllar gas 
- see Sect. 4 and Paper II). Therefore we find no strong 
evidence for binarity. 

5.3. Radio and X-ray emission properties 

36 Lyn was observed at 6 cm using the VLA by Linsky 
et al. (1992), who found an upper limit (3a) of 0.62 mJy. 
More recently, it has been detected (at phase 0.73 accord- 
ing to our ephemeris) with a flux of 0.45 ± 0.05 mJy at 3.6 
cm (S. Drake, private communication). 

Berghofer et al. (1996) found that ROSAT obtained no 
detection of X-rays at the position of 36 Lyn, reporting 
logT, < 29.84. 

5.4. Magnetic field geometry 
5.4.1. Longitudinal field variation 

For a tilted, centred dipole of polar strength B p , the varia- 
tion of the longitudinal magnetic field (B z ) with rotational 
phase 4> is given by Preston's (1967) well-known relation 

15 ± u 

(B z ) = Bp-^^ — ^y(cos/3cosi±sin/3sinicos27r(0-(/)o)), (2) 

where u denotes the limb darkening parameter and 0o the 
phase of longitudinal field maximum (equal to u = 0.4 and 
approximately <fio = 0.25 in the particular case of 36 Lyn). 
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Fig. 6. A comparison between the observed flux distribution in 36 Lyn (filled squares) and the theoretical distributions 
calculated with three selected atmospheric models: 13000g40pl0k2 (dotted line), 13300g40p05k2 (thick line, and the 
adopted model), and 13600g37p00k2 (thin line). A comparison between the observed H/3 (left panel) and Ha (right 
panel) line profiles and synthesized with f 3300g40p05k2 (full thick line) and f 3600g37p00k2 (dotted line) models. 



The inclination and obliquity angles i and ft are related 

by 



n 1 ~ r 

tan p — cot i, 

I + r 



(3) 



where r = —0.83 ± 0.1 is the ratio of the longitudinal field 
extrema according to the best- fit sine curve . 

From r ~ -1.0 we obtain straightforwardly that either 
i or (3 is close to 90°, a result consistent with the rigid 
rotation model developed considering the stellar radius, 
period and usini. 

For our best-fit inclination i — 90°, (3 cannot be deter- 
mined from Eq. (3) and Eq. (2) reduces to 



(B z y 



= B r 



15 



u 



20(3 - u ) 



sin/3. 



(4) 



4 Fig. 3 suggests that there may exist a systematic difference, 
at the level of a few hundred G, of the longitudinal field vari- 
ation near magnetic minimum as diagnosed using metal line 
versis Balmer line measurements. Such differences frequently 
exist, are poorly understood, and may be related to the pres- 
ence of nonuniform surface abundance distributions. In the case 
of 36 Lyn, sine fits to the entire dataset or to only those mea- 
surements obtained from H/3 provide r values that are identical 
to within the uncertainties. 



In other words, the magnetic obliquity /3 and polar 
strength of the dipole B p cannot be disentangled when i = 
90° , and we can only determine their product B p sin j3 = 
3380 ± 170 G. 

In the case of our minimum acceptable inclination i = 
56°, we obtain using Eq. (3) that (3 = 82 ± 5°, and Eq. 
(2) gives B p = 3760 ± 170 G. Therefore, considering the 
full range of inclinations 56° — 90° admitted by the data, 
we infer from the longitudinal field variation that 3210 < 
B p < 3930 G. 



5.4.2. Stokes V profiles 

Using the method described by Donati et al. (2001) and 
Wade et al. (2005), we have compared the phase variation 
of the Stokes V profiles with the Zeeman2 (Landstreet 
1988, Wade et al. 2001) predictions of a grid of 9000 dif- 
ferent magnetic field models, obtained by systematically 
varying i, Bd and (3 in the range 56-90°, to -4000 G, 
and to 180°, respectively. For each model, the reduced 
X 2 of the observed and calculated Stokes V profiles was 
calculated. 
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The best-fit solutions for the dipolar strength B d and 
the magnetic obliquity /3, for all admissible inclinations, 
are B d = -1300 G and /3 = 100. Although both methods 
recover effectively the same value of the obliquity (/? = 80° 
and (3 = 100° produce identical longitudinal field varia- 
tions), the dipole strength obtained from the profile fits 
is abour 3 times smaller than that determined from the 
longitudinal field. 

We strongly suspect that these (significant) differences 
result from the important modulation of the Stokes V pro- 
files by the non-uniform distribution of Fe on the stellar 
surface. In fact, the best-fit model for Stokes V, while re- 
producing the amplitudes of the weak signatures around 
phases 0.9-0.2 rather well, is totally unable to reproduce 
the stronger Stokes signatures at phases 0.5-0.8. This 
failure is reflected in the very high reduced y 2 statis- 
tic {x 2 / v ~ 10.4) obtained for the Stokes V profiles. 
Moreover, the structure and radial velocities of features 
apparent in the observed Stokes V profiles (see Fig. 1) 
are not evident in the calculated profiles. This suggests 
that the model we have used to calculate the profiles of 
the LSD mean profiles is an obvious oversimplification of 
the true surface structure of 36 Lyn. Additional attempts 
to reproduce this structure using higher-order multipolar 
magnetic field models were unsuccessful. We therefore con- 
clude that the structure of the Stokes V profiles is likely 
due to a strongly nonuniform distribution of Fe on the 
surface of the star. This will be discussed further in Sect. 
6. 

Finally, we report that Stokes Q and U profile cor- 
responding to the best-fit Stokes V dipole model (with 
| .Bp | = 1300 G) have amplitudes similar to the typical 
noise in the associated LSD profiles. For the dipole model 
derived from the longitudinal field data (with \B p \ = 
3380 G), the predicted linear polarisation amplitudes are 
several times larger than the observed Stokes Q and U 
upper limits. This is consistent with results obtained for 
other magnetic Ap/Bp stars by Wade et al. (2000a), and 
may reflect the influence of non-uniform abundance dis- 
tributions, or small-scale structure of the magnetic field, 
on the disc-integrated polarimetric signal. 

6. Photospheric chemical abundances 

6.1. Mean abundances 

Mean (phase-averaged and surface-averaged) photospheric 
abundances were determined from analysis of both optical 
and UV spectra of 36 Lyn. Some previously published UV 
abundances were also obtained from the literature. All 
abundances are reported in Table 5. 

6.1.1. Optical abundances 

Synthetic spectrum calculations for the whole spectral re- 
gion covered by the optical observations (4500-6600 A) 
were performed with the synthetic spectrum codes Synth 
(Piskunov 1992) and STARSP (Tsymbal 1996). Selected re- 



gions of the spectrum were also modeled using Zeeman2 
(Landstreet 1988, Wade et al. 2001). For hydrogen line 
profiles we used the more recent calculations of the hy- 
drogen line opacities (Barklem at al. 2002). We did not 
include the magnetic field for SYNTH or STARSP synthetic 
calculations; instead, we employed a pseudomicroturbu- 
lence of 2 kms -1 to account for magnetic intensification. 
The influence of a dipolar magnetic field of surface inten- 
sity 3 kG was employed for the Zeeman2 calculations. 
The adopted model atmosphere is an ATLAS9 model 
with the adopted parameters described in Sect. 5. We 
have explored the range of model parameters discussed in 
Sect. 5.1, derived from the photometry, energy distribu- 
tion and Balmcr lines, and we find that the detailed choice 
of model atmosphere has a negligible impact on the de- 
rived abundances. With a few exceptions, all atomic data 
were extracted from the VALD database (Kupka et al. 
1999, Ryabchikova et al. 1999). Rare earth abundances (Pr 
and Nd) were derived from the lines of the second ioniza- 
tion stage with the transition probabilities from Cowley & 
Bord (1998), Bord (2000) - Nd III, and from Bord (private 
communication) - Pr III. Two lines of Ga II, AA 6334, 6419 
were used for the gallium abundance determination in the 
optical region. Hyperfine and isotopic splitting (Karlsson 
& Litzen 2000) were taken into account. 

In order to fit the line profiles, it was necessary to 
determine the projected rotational velocity v sini. Various 
estimates of wsinz exist in the literature, ranging from 
40 kms _1 (Abt & Morrel 1995) to 60 kms _1 (Abt, Levato 
& Grosso 2002). We find vsmi = 48 ±3 kms -1 , consistent 
with the recent results of Royer et al. (2002) who find 
■usini = 49 kms -1 . 

Fig. 8 shows a part of the 36 Lyn optical spectrum and 
synthetic spectrum fits. 

6.1.2. Ultraviolet abundances 

Ultraviolet abundances were inferred from spectrum syn- 
thesis using the SYNSPEC code (Hubeny, Lanz & Jeffery 
1994), assuming a Kurucz (1993) ATLAS9 model atmo- 
sphere. 

Notably, we discovered that the silicon abundances de- 
rived from Si II and Si III lines differed, a problem that we 
ultimately traced to a series of blends in the Si in lines 
and which may potentially lead to erroneous abundances 
of Si derived for Bp stars from IUE data (two good, un- 
blended Si in lines are those at 1370.0 Aand 1417.2 A.) We 
were finally able to distinguish the unblended and blended 
Si in lines by their variability, as the former show varia- 
tions in their radial velocity with phase, as reported from 
an optical Si line study by Stateva (1997). 

6.2. Mean abundances and evidence for vertical 
abundance non-uniformities 

In our analysis of optical spectra of 36 Lyn, we find no 
systematic differences in the abundances required to fit 
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Fig. 8. A comparison between the observed and synthesized spectra in the regions of Alii, Con, and Nin lines (left 
panel) and Sin, Gail, and Ndin lines (right panel). Calculations with the Al, Co, Ni, Ga and Nd abundances derived 
from our analysis are shown by the full line, and those with the UV abundances (except Fe and Ti) are shown by the 
dotted line. The dotted line also represents the Ndlll line with a 1.0 dex reduced abundance. 



36 Lyn Optical 1 36 Lyn UV 
Element logn/n to t logn/ntot logn/ntot 

-1.05 
-3.48 
-3.11 
-3.95 
-5.71 
-4.46 
-5.57 
-4.49 

-6.59 
-4.83 
-6.54 
-7.05 
-6.37 

-6.65 
-4.37 

-7.12 
-5.79 
-7.83 
-9.16 
11.33 
10.54 



Table 5. Mean photospheric chemical abundances of 36 
Lyn. Also shown are the respective solar abundances. 
Sources: ^his work, 2 Smith (1993), 3 Smith & Dworetsky 
(1993), 4 Smith (1994), 5 Smith (1996). 



strong versus weak lines of elements with sufficiently rich 
line spectra. This suggests that stratification of these ele- 
ments in the atmosphere of 36 Lyn is either very weak or 
nonexistent. 

For a few elements, the abundances obtained from op- 
tical lines differ from those derived from UV lines. These 



differences cannot be explained by small differences in the 
adopted atmospheric parameters. The largest disagree- 
ments occur for Mn, Co, Ni and Ga which appear to be less 
abundant from the UV analysis, while Al is opposite. The 
Ga problem is well known (Lanz et al. 1993). Dworetsky & 
Smith (1993) derived abundances from the low excitation 
UV Mn II, Co II and Ni n lines, while our analysis is based 
on high excitation lines. Therefore a difference in the UV 
and optical abundances may reflect stratification of these 
elements, although the discrepancy could be attributed to 
unidentified blends given the relative scarcity of lines of 
these elements in the optical. An examination of the oscil- 
lator strengths used in the various studies indicates that 
systematic errors in gf values cannot explain the differ- 
ences. For Co II lines, neglect of hyperfine structure may 
be an important contributor. 

Overall, we find that He is depleted relative to solar, 
while Fe peak elements are strongly enhanced, by one dex 
or more. Abundances of two rare earth elements (REEs, 
Pr and Nd) are found to be enhanced by approximately 3 
dex. Mn, Co, Ni and Ga may be stratified. 

6.3. Horizontal abundance non-uniformities 

The variations of lines of a number of elements in the high- 
resolution optical spectrum of 36 Lyn as well as the LSD 
profiles are suggestive of non-uniform horizontal distribu- 
tions of their abundances. However, some care is required 
in interpreting these variations: Smith & Groote (2001) 
caution that the magnetically-confined winds of Bp stars 
can produce optical line profile variability due to occulta- 
tion of the stellar disc by the circumstellar material, and 
that in a number of hot Bp stars a significant component 
of the optical metallic-line variability can be attributed to 
this phenomenon. 

In the spectrum of 36 Lyn, weak variability occurs in 
lines of He, Fe, Si, Cr, Ti, and Mg (based on our own 



He -1.65 ±0.05 

C -3.0 - -2.7 -3.48 1 

O -3.4 - -3.0 

Ne -3.8 - -3.95 

Na -5.30 

Mg -4.35 -4.44±0.22 2 

Al -6.30 -6.10±0.15 2 

Si -4.10 ±0.2 -4.5 : ±0.4 X 
-4.70±0.10 2 

P -6.30 

S -4.8 - -5.0 

CI -5.2 

Ti -5.00 ±0.25 

Cr -5.00 ±0.3 -5.3 1 , 

-5.3±0.4 3 

Mn -5.8 -6.80±0.15 3 

Fe -3.45 ±0.25 -3.2 1 
-3.90 ±0.1 3 

Co -5.00 ±0.15 -6.2±0.2 3 

Ni -4.90 ±0.15 -5.7±0.1 3 

Cu -8.1±0.6 4 

Ga -6.05 ±0.1 -6.90±0.55 5 

Pr -8.1 ±0.2 

Nd -7.7 ±0.2 
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Fig. 9. Fe Doppler image from Fe II A4522. Right frame - comparison of observed and calculated profiles. Left frame 
- Derived abundance map in rectangular projection. Tick marks indicate phases of 36 Lyn observations used to derive 
the map. 



spectra, as well as the results of Stateva 1997). Some of 
these elements show extrcma of their equivalent width 
phase variations which coincide with the magnetic ex- 
trema (phases 0.25/0.75, either in phase or anti-phase, e.g. 
He, Fe) whereas others show variations with extrema at 
crossover phases (phases 0.0/0.5, e.g. Si). We have also ex- 
amined the dependence of EW variation shape and phase 
of Fe lines as a function of excitation potential, and find no 
significant dependence. Based on (1) the variety of equiv- 
alent width variability characteristics, (2) the large range 
in excitation potential of lines that undergo variability, (3) 
the modulation of the shapes of Stokes V profiles and (4) 
the apparent difference in the extrema of the longitudinal 
field variations as diagnosed in metal versus H lines, we 
conclude that the optical line profile variations result sub- 
stantially from surface structure (magnetic + abundance 
non- uniformity) , rather than being due to the circumstel- 
lar material. This conclusion is supported by an examina- 
tion of the high-resolution line profiles. As can be seen in 
Figs. 1 and 9, Fe lines exhibit variations in Stokes I and 
V profile shapes, consistent with Doppler modulation by 
surface non-uniformities. 

Stateva (1997) attempted an harmonic reconstruction 
of the surface distributions of He, Si, Cr and Mg using 
equivalent width variations extracted from photographic 
spectra of 36 Lyn, and assuming that the variations were 
due entirely to surface non-uniformities. She was able to 



reproduce the (relatively low signal-to-noise ratio) equiva- 
lent width and radial velocity variations, and the resulting 
maps suggest that He, Si and Cr are concentrated primar- 
ily in the plane of the magnetic equator. 

Here, we apply the more stringent test of attempting 
to reproduce the detailed line profile variations using the 
Doppler Imaging technique, for the particular case of Fe. 
Because we have only 11 high-resolution spectra, we are 
not attempting to derive a definitive map of the surface 
Fe distribution. Our primary goal is rather to demonstrate 
that the detailed profile variations of moderate and high- 
excitation Fe lines in the optical can be substantially ex- 
plained by an inhomogencous surface distribution of this 
element, and are not significantly affected by absorption 
due to the circumstellar material. 

Doppler imaging is a technique which allows the in- 
version of observed spectral line profile variations to pro- 
duce a two-dimensional map of the abundance distribu- 
tion at the surface of a star. An inhomogeneity on the 
stellar surface leads to distortions in the spectral line pro- 
files. These distortions appear on the blue side of the line 
profile as the spot comes into sight on the visible hemi- 
sphere of the star, travel across the spectral line and dis- 
appear on the red wing as the spot rotates out of sight. 
The longitudinal position of the inhomogeneity is obtained 
directly from the velocity position of the distortion within 
the profile, whereas its latitude has to be deduced from 
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spectral time series. In INVERS12, the Doppler Imaging 
code we applied, developed by N. Piskunov and refined by 
O. Kochukhov (Kochukhov et al. 2004b), the equation of 
radiative transfer is solved in a way that allows the calcu- 
lation of maps of several chemical elements simultaneously 
and self-consistently. In this program the precalculated lo- 
cal line profile tables used in former versions of the code 
are replaced by the calculation of specific intensities for 
each visible surface element on each iteration. 

For 36 Lyn, we assumed various inclination angles be- 
tween 56° and 90° to perform the Doppler Imaging. The 
best fit to the profiles was found for i — 60°, although 
the dependence of the fit quality on i was rather weak 
(and in particular this dependence cannot be used to con- 
strain the inclination). We derived the Fe abundance on 
the surface of the star using several Fe II lines, including 
A4520, A4522, A5018, A5534 and A5780, each of which pro- 
duced a similar Fe distribution. In Fig. 9 we show the Fe 
distribution, in the rotational frame of reference, recon- 
structed from Fen A4522.6 (% e = 2.8 eV). The distribu- 
tion is characterised by a very patchy ring-like structure 
between -30°and +30° in stellar latitude, with a total con- 
trast of somewhat more than 3 dex. Due to the uncertainty 
in the magnetic field geometry, it is difficult to conclude 
whether any obvious geometrical relationship exists be- 
tween the Fe surface distribution and the magnetic field. 
Our success in interpreting the profile variations of the 
Fe II lines supports our interpretation of the Fe (and other 
metal) line variations are due to Doppler modulation of 
surface structure. The details of the derived map should 
be interpreted conservatively, given that the map was de- 
rived using a single line with relatively coarse rotational 
phase coverage. 

7. Discussion and conclusion 

We have performed a thorough reevaluation and analy- 
sis of the physical and spectroscopic characteristics of the 
B8p He-weak star 36 Lyncis. Based on measurements of 
the longitudinal magnetic field, as well as measurements 
of the core equivalent width of the Ha line, we have rede- 
termined the rotational period P ro t = 3.83475 ± 0.00003 
days. An examination of the variability of the Ha core 
equivalent width reveals strong changes coincident with 
magnetic crossover phases, which we interpret as due to 
occultation of the stellar disc by magnetically-confined cir- 
cumstellar gas. We have determined the effective temper- 
ature, radius, mass, age, projected and true rotational ve- 
locities, the rotation axis inclination, and we have sum- 
marised evidence for binarity, radio and X-ray emission. 

We have found systematic and significant discrep- 
ancies between the temperature and gravity determined 
from the optical and UV energy distribution and from 
the Balmer lines, concluding that they probably show 
true limitations of the ATLAS9 model atmospheres in de- 
scribing the real photospheric structure of 36 Lyn. The 
ATLAS9 models we have employed are calculated using 
scaled solar metallicities, while individual abundances in 



the atmosphere of 36 Lyn show large diversity. To check 
the influence of anomalous abundances on the model at- 
mosphere structure, a few model atmospheres were calcu- 
lated with LLmodels model atmosphere code (Shulyak 
et al. I2004|l which takes into account the contribution of 
individual line opacities to the total continuum opacity 
through direct synthetic calculations. This code allows us 
to examine the influence of the particular abundance table 
of 36 Lyn on the model atmosphere structure. The flux 
distribution for an LLmodels model with T cff =13300, 
log g = 4.0 and the abundances of 36 Lyn actually provides 
a somewhat better simultaneous fit to both the energy 
distribution and Balmer lines. We conclude that a model 
atmosphere with individualised abundances improves the 
overall fit to the observables, although discrepancies are 
still evident (a surface gravity below 4.0 is still required 
to fit the H lines) . Possibly we are seeing the influence of 
peculiar and possibly stratified abundances on the photo- 
spheric temperature and pressure structure - effects which 
are not included in LLmodels. 

We have modeled Least-Squares Deconvolved Stokes 
V profiles and the longitudinal field variation to constrain 
the magnetic field geometry. We have analysed both op- 
tical and UV spectra and, using spectrum synthesis tech- 
niques, we have determined the abundances of 21 ele- 
ments. We find no evidence for vertical stratification of 
most elements, although for a few elements (e.g. Mn, Co, 
Ni, Ga) we have derived abundances in the optical which 
differ from those in the UV, suggesting that these elements 
may be stratified. The influence of using an LLmodels 
individualised abundance model atmosphere on inferred 
abundances is small (±0.1 dex) for all elements except He, 
for which the abundance is decreased by 0.4 dex. On the 
other hand, using the individualised abundance model de- 
creases by 0.15 dex the difference between the optical and 
UV Ga abundances (Ga optical abundance is decreased 
by 0.1 dex and Ga UV abundance increased by 0.05 dex). 

We find convincing evidence exists for nonuniform sur- 
face (horizontal) distributions of some elements. Using the 
Doppler Imaging technique, we have constructed an ap- 
proximate surface map of the Fe abundance distribution, 
which is characterised by a patchy ring-like structure be- 
tween —30° and +30°, and a contrast of somewhat more 
than 3 dex. The results reported in this paper set the stage 
for a detailed study of the structure and geometry of the 
circumstellar material in Paper II. 
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